CROSS DIFFUSION AND NONLINEAR DIFFUSION PREVENTING 
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Abstract. A parabolic-parabolic (Patlak-) Keller-Segel model in up to three space di- 
mensions with nonlinear cell diffusion and an additional nonlinear cross-diffusion term is 
analyzed. The main feature of this model is that there exists a new entropy functional, 
yielding gradient estimates for the cell density and chemical concentration. For arbitrar- 
ily small cross-diffusion coefficients and for suitable exponents of the nonlinear diffusion 
terms, the global-in-time existence of weak solutions is proved, thus preventing finite-time 
blow up of the cell density. The global existence result also holds for linear and fast diffu- 
sion of the cell density in a certain parameter range in three dimensions. Furthermore, we 
show bounds for the solutions to the parabolic-elliptic system. Sufficient conditions 
leading to the asymptotic stability of the constant steady state are given for a particular 
choice of the nonlinear diffusion exponents. Numerical experiments in two and three space 
dimensions illustrate the theoretical results. 



1. Introduction 

Patlak [26] and Keller and Segel [19] have proposed a partial differential equation model, 
which describes the movement of cells in response to a chemical signal. The cells move 
towards regions of higher signal concentrations. As the cells produce the signal substance, 
the movement may lead to an aggregation of cells. The more cells are aggregated, the more 
the attracting chemical signal is produced by the cells. This process is counter-balanced by 
cell diffusion, but if the cell density is sufficiently large, the nonlocal chemical interaction 
dominates and results - in two and three space dimensions - in a blow up of the cell 
density (see the reviews [12, 15] for details). Denoting by p = p{x,t) the cell density and 
by c = c{x, t) the concentration of the chemical signal, the Keller-Segel model, in its general 
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form, is given by 

dtp = dw{D{p)Vp - x(p)Vc) + c), 
adtc = Ac + i?2(p, 5*), X eVt, t> 0, 

where C M"^ (c? > 1) is a bounded domain, D{p) is the cell diffusivity, x{p) the chemotac- 
tic sensitivity, and -Ri(p, c) and -R2(p, c) describe the production and degradation of the cell 
density and chemical substance, respectively. Here, a = corresponds to the parabolic- 
elliptic case and a = 1 to the fully parabolic problem. The equations are supplemented by 
homogeneous Neumann boundary and initial conditions: 

D(p)(Vp ■ z/) = Vc- z/ = on t > 0, 
p(-,0)=po, ac(-,0) = aco in il, 

where v denotes the exterior unit normal to the boundary dVt (which is assumed to exist). 
The classical Keller-Segel model consists in the choice -D(p) = 1, x(p) = p, Ri{p,c) = 0, 
and R2{p, c) = p — c. 

Motivated by numerical and modeling issues, the question how blow up of cells can 
be avoided has been investigated intensively the last years. Up to our knowledge, four 
methods have been proposed. In the following, we review these methods. 

The first idea is to modify the chemotactic sensitivity. Supposing that aggregation stops 
when the cell density reaches the maximal value poo = 1, one may choose x{p) = 1 ~ P- 
this volume-filling case, the cell density is bounded, < p < 1, and the global existence 
of solutions can be proved [9]. Furthermore, if x(p) = P^ with < (3 < 2/d, the solutions 
are global and bounded, thus preventing finite-time blow up [16]. Global solutions are also 
obtained when the sensitivity depends on the chemical concentrations in an appropriate 
way, see, e.g., [1, 13]. 

A second method consists in modifying the cell diffusion. In the context of the volume- 
filling effect. Burger et al. [3] suggested the cell equation dtp = div(p(l — p)V(p — c)). Then 
the parabolic-elliptic model possesses global solutions. Global existence results can be 
achieved by employing the nonlinear diffusion D{p) = p", which models the local repulsion 
of cells. When fl'{D{s)/s)ds grows faster than logp for large p, a priori estimates showing 
that solutions are global and uniformly bounded in time were obtained in [4, 21]. Adding 
the nonlinear sensitivity x{p) = with a > 1 and 2 < /3 < a + 2/d, global existence 
results were achieved in [17]. The solutions are uniformly bounded in time if a > 2 — 4/ d 
[22]. The existence of global bounded classical solutions to a fast- diffusion Keller-Segel 
model with -D(p) = (1 — p)~", where a > 2, has been proved in [6]. The same result holds 
true when we choose x(p) = p(l — pY with /3 > 1 — a;/2, and the solution is still global in 
time (but possibly not classical) if /3 > 1 — a [30] . 

A third approach is to consider nonvanishing growth-death models Ri ^ 0, since one 
may expect that a suitable death term avoids cell aggregation. Indeed, taking i?i(p, c) = 
p(l — p)(p — a) for some < a < 1, the global existence of solutions is proved in [2]. 
In the logistic-growth model i?i(p,c) = p(l — p^'^), a global weak solution exists for all 
7 > 2 — l/c? [29]. These results have been obtained for the parabolic-elliptic model. 
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Recently, a fourth way to obtain global existence of solutions has been proposed [14]. 
The idea is to add a cross- diffusion term in the equation for the chemical signal: 

dtp = div(Vp - pVc), 
adfC = Ac + 6Ap + p — c in fi, t > 0, 

where 6 > 0. At first sight, the additional cross-diffusion term seems to cause several 
mathematical difficulties since the diffusion matrix of the above system is neither symmetric 
nor positive definite, and we cannot apply the maximum principle to the equation for the 
chemical signal anymore. All these difficulties can be resolved by the observation that the 
above system possesses a logarithmic entropy, 

f I c^l 
Eo{p,c)= p(logp-l) + a— dx, 

t/ _ 

allowing for global existence results and revealing some interesting structural properties of 
the system. In fact, the entropy production equation 

and suitable Gagliardo-Nirenberg estimates for the right-hand side lead to gradient es- 
timates for and c. Another motivation for the introduction of the additional cross 
diffusion is that, whereas finite-element discretizations of the classical Keller-Segel model 
break down some time before the blow up, the numerical solutions to the augmented 
model exists for all time, which may lead to estimates of the blow-up time. This question 
is currently under investigation. 

In [14], the existence of global weak solutions has been proved in the two-dimensional 
situation only. In this paper, we generalize this result to three space dimensions by allowing 
for nonlinearities in the cell diffusion terms. Since nonlinear diffusion in the cell equation 
helps to achieve global existence results (see above), we suggest, in contrast to [14], a 
nonlinear cross-diffusion term. More precisely, we consider the equations 

(1) 9tp = div(V(p"^)-pVc), 

(2) adtc = Ac + 5A(p") + p - c in fi, t > 0, 
subject to the no-flux and initial conditions 

(3) (V(p™) - pVc) ■ z/ = V(c + 5p") ■ z/ = ondn, t>0, 

(4) p(.,0) = po, ac(-,0) = aco in fi. 

Notice that these boundary conditions are equivalent to Vp ■ u = Vc ■ z/ = on dQ for 
smooth positive solutions. 

In two space dimensions, the case m = n = 1 is covered by [14]. If m > 3 — 4/d, 
2 < d < 3, the nonlinear diffusion already prevents blow-up of the solutions without 
additional cross diffusion, see [20, 21, 22]. The question remains if we can allow for linear 
and fast diffusion of cells, m < 1, for some n > 1, and still obtain global existence results. 
In this paper, we show that this is indeed true. For instance, we show that in the presence 
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of the additional cross diffusion term and in three space dimensions, we can allow for the 
classical cell diffusion exponent m = 1 and still obtain global existence results. This shows 
that the result of [14] can be generalized to the three-dimensional case if the cross diffusion 
is of degenerate type. These remarks motivate us to restrict ourselves to the case m > 
and n > 1. Our first main result is as follows. 

Theorem 1 (Global existence). Let C M"' (1 < (i < 3) be a bounded domain with dQ G 
C^'^. Let a > 0, m > 0, n > 1, and let p = {m + n — l)/2 satisfy 1 — n/d < p < min{m,n}. 
Furthermore, let < po G and aco G Then there exists a global weak solution 

{p,c) to (l)-(4) satisfying p>0 in Q, t > 0, and, for some s G (1,2], 



where Q = n/d + p > 1. 

Remark 2. A weak solution is to be understood in the standard sense by testing the system 
of equations against compactly supported smooth functions in C°°((0,T) x Q)). Due to 
the regularity properties of the solution, however, test functions in L'*(0, T; iy^'*(r2)) are 
sufficient for the weak formulation of both equations in the fully parabolic system to be 
well defined. For the parabolic-elliptic system, we show in Section 4 that we can even allow 
for test functions in L^{0,T] H^{Q)). □ 

Let us discuss the conditions on p which are equivalent to 

2 

(5) m — l<n<m + l, m + n + -n > 3. 

The areas of admissible values for {m,n) are illustrated in Figure 1. Notice that the bands 
between n — l<m<n + l continue to the right. 

n d=l n (1,2 n d=3 



p G oo; L"m n 4^(0, oo; L^^m, 

p'", p" G LUO, oo; W''%n)), pVc G ^^,(0, oo; L%n)), 

ac G (0, oo; I'm, c G LL(0, oo; H\n)), 

dtp, adtceLlMoo;{W'''m'), 
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diffusion, m = 1, the above conditions are satisfied for any 1 < < 2 (if = 2) and 
6/5 < n < 2 (if d = 3). Hence, tlie degenerate cross- diffusion term prevents blow up in 
finite time even in tlie case of linear cell diffusion in three dimensions. In short, one of the 
conditions in (5) is needed to derive a nice bound on an entropy functional and the others 
for suitable compactness and continuity properties of the approximated sequences. 

The key idea of the proof of Theorem 1 is the observation that system (l)-(2) possesses 
an entropy functional, 

(6) E{p,c)= [ (-^ + a^]dx, 



.n-1 26. 

useful to derive a priori estimates. Indeed, differentiating formally this functional, we 
obtain the entropy production equation 

dE f fmn._ 1 , 1 



dt 



recalling that p = (m+n — 1)/2. We will show in the proof of Lemma 10 that the right-hand 
side can be estimated for any /3 > as follows: 

(7) / pcdx <P [ \Vff\'dx + C{/3, \\p\\mn)) + \ ! {\Vc\^ + ^)dx, 

Jn Jn ^ Jn 

under the restriction 1 — 2/(i < p, which follows from the conditions in (5) for 1 < c? < 3. 
The assumptions m > and n > 1 imply that p > 0, thus only for c? = 3 we obtain the 
restriction p > 1/3. Let us remark that in the case c? = 3, the restriction 1 — 2/c/ < j9 is 
redundant, that is, conditions (5) together with n > 1 and m > imply that p + 2/d > 
p + n/ d > 1 for 72 G (1, 2], as well as p > 1 since 2p = m + n — 1 >2for?2>2 with 
m > n — 1 > 1. 

The existence of the entropy functional (6) implies the existence of so-called entropy vari- 
ables which makes the new diffusion matrix positive (semi-) definite. Indeed, introducing 
the entropy variables 

dE n dE c 



dp n — 1^ ' dc (5' 

system (l)-(2) can be written as 

In hyperbolic or parabolic systems, the existence of an entropy functional is equivalent 
to the existence of a change of unknowns which "symmetrizes" the system [8, 18]. (For 
parabolic systems, "symmetrization" means that the transformed diffusion matrix is sym- 
metric and positive definite.) In system (8), the diffusion matrix is nonsymmetric, but still 
positive semi-definite. 

The existence proof is based on the construction of a problem which approximates (8). 
First, we replace the time derivative by an implicit Euler approximation with time step 
r > and add a weak form of the fourth-order operator ^(A^r — div(| Vr pVr) + r) {e > 0) 
to the first component of (8), which guarantees the coercivity of the elliptic system in 
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H'^(Q) with respect to r. The existence of weak approximating solutions (r^jb^) is shown 
by using the Leray-Schauder fixed-point theorem. At this point, we need the restriction 
p < m, which is equivalent to n < m + 1, in (5) to ensure the continuity and coercivity. 
The discrete entropy estimate implies a priori estimates uniform in the approximation 
parameters r and e, which allow us to pass to the limit (r, e) — )■ 0. 

There are two technical difficulties in the limiting procedure. The first one is that the 
entropy estimate yields a uniform bound for pP in H^{Q), but an estimate for dtPs in 
{H^{Q)y. If p < 1, this implies a bound for in W^''^{Q) for some r > 0, and we can 
apply the Aubin lemma to conclude the relative compactness of the family (pe)e>o- If p > 1, 
we infer this property using a variant of the Dubinskii lemma (see Lemma 7). The second 
difficulty is to ensure the strong convergence of the family (p£)e>o in L'^{^ x (0,T)). In 
two space dimensions d = 2 (and with n = m = 1), this has been proved in [14]. However, 
for (i = 3 (and p < 1), our uniform estimates in Lemma 12 need additional assumptions 
on the diffusion parameters, namely p > 1 — n/d and p < n, or equivalently, the remaining 
two conditions in (5): m + n + 2n/d> ?> and m — 1 < n. 

Our second main result concerns some qualitative properties of the solutions to (l)-(4) 
using the entropy functional. First, we prove L°° bounds for the solutions to the parabolic- 
elliptic system generalizing the results of [21, 22] to this situation. 

Theorem 3 (Boundedness in Let the assumptions of Theorem 1 hold and let a = 0. 

Then, for any T > 0, the solution {p,c) to the parabolic- elliptic system (l)-(4) satisfies 

||p||L°°(0,T;L°°(n)) + ||c||L=o(0,T;L°°(n)) < C (T) , 

where the constant C{T) > depends on T > 0. 

Second, we are able to show the asymptotic stability of solutions to the constant steady 
state. Due to the special structure of the entropy functional, we can allow for a very 
particular choice of the parameters m, n, and 6 only. 

Proposition 4 (Long-time decay for m = 1, n = 2). Let C M'^ (1 < d < 3) be a bounded 
domain with dQ G C^'^. Let po G L°°{Q), m = 1, n = 2, and 6 > Cp/A, where Cp is the 
constant of the Poincare inequality in L'^{Q). Then the solution to the parabolic- elliptic 
system (32)-(33) with a = 0, constructed in Theorem 15, decays exponentially fast to the 
homogeneous steady state in the sense that 

\\p{;t)-p*U2^n) < Ce-^\ ||c(-, t) - c* < Ce-^\ 

where C > is some constant and k, = min{l,45 — C|,}/(45). Moreover, any smooth 
solution {p,c) to the fully parabolic system (32)-(33) with a = 1 has the decay properties 

M;t)-p*\\L2(^n) < Ce-^\ \\c{;t)-c*\\LHn) < Ce"'^* 

for all t > 0. 

The paper is organized as follows. In Section 2, we prove an inequality which is needed for 
the proof of (7) and we show a compactness result which combines the lemmas of Aubin 
and Dubinskii. Theorems 1 and 3 are shown in Sections 3 and 4, respectively, whereas 
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Proposition 4 is proved in Section 5. Finally, in Section 6 we present some numerical 
results in two and three space dimensions which illustrate the effect of the exponent n. 

2. Auxiliary results 

Lemma 5. Let Vt dW^ [d > 1) he a hounded domain with dVt G Furthermore, let 

f3 > and 

• ifd> 3; either I - 2/d <p <l and q = 2d/{d + 2) 

or p > 1 and q = p + 1/2, 

• if d < 2: either < p < 1, q > 1, and p+l/q > 3/2 - 1/d 

or p > 1 and q = p + 1/2. 

Then there exists a constant C{(3, ||p||Li(n)) > such that for all p G L^{Q) satisfying p > 
in Q and fP G H^{VL), the following inequality holds: 

\\p\\Un)<mf^\\hin) + C{PA\p\\LHn)). 

Notice that the continuous embedding H^{Q) ^ L^{Q), where 1 < s < 2d/{d — 2) if 
d > 3 and l<s<cx)if(i<2, shows that pP G H^{il) implies that p G L'^P(f2), and the 
condition q < sp has to be imposed. This condition is satisfied for the above choices of p 
and q. 

Proof. First, let < p < 1 and (p, q) be given as in the Lemma. The Gagliardo-Nirenberg 
inequality, see e.g. [10, Theorem 10.1] and [31, Theorem 1.1.4], gives 

where 6 = dp{l — l/g)/(l — d/2 + dp) and (7 > is here and in the following a generic 
constant. The conditions p>l — 2/(iif(i>3 and g > 1 if ci < 2 imply that 6 > 0. For 
all space dimensions, it holds that p + 1/q > 3/2 — 1/d which is equivalent to 6' < p < 1. 
Then the inequality 6/p < 1 allows us to apply the Young inequality: 

\\p\\Un)<mp'\\Un) + Cm\p\\l^^^^^^^^ 

proving the first case. 

Next, let p > 1 and q = p + 1/2. Notice that the Poincare inequality implies that 



L2(Q) 



< 



f- fdx 



n 



+ \\f\\LHn)<Cp\\Vfh^n) + \\f\\LHn). 



This together with the Holder inequality leads to 

iipiiL(n) = m%u = ii/p^/iK^(.) 

< \\Pl%\nM^rAn, < ^(I|V/||K^(.) + ll/llK'(n))IIP^/^ilK^(.) 
= C\\'^P^\\L2\n)\\p\\i!i\n) + C'IIPllLp'(f2)IIPllL'i^(n)- 



8 J. A. CARRILLO, S. HITTMEIR, AND A. JUNGEL 

Furthermore, using the interpolation inequahty with 1/p = 6/q + {l — 6)/l or, equivalently, 
p9/q = {p-l)/{q-l)>0, 

WPWiiin) ^tv||vp \\L2(^n)\\P\\L^n) + '-'\\P\\L1{n)\\P\\L^n) \\P\\L^n)■ 
Since g > 1, we may employ the Young inequality for the first summand to obtain 

\\p\\Un) < /3||V//||i.(.) + C(/3)||p||i/(-/) + C||p||if/^^^ 

Since 1 < p < g, it follows that 2p6/q < 2, which allows us to use the Young inequality for 
the second summand: 

MUn) < mPlL^n) + Cm\p\\LHn''' + l\\p\\Un) + CW, Mmn))- 

The lemma is proved. □ 

Next, we recall a compactness result. Let {ahp){x,t) = p{x,t — h), t > h > 0, he a shift 
operator. 

Lemma 6 (Dubinskii). Let Q G M.'^ {d > 1) be a bounded domain with dQ G C^'^ and let 
T > 0. Furthermore, let p > 1, q > 1, and s >0, and let (p^) be a sequence of nonnegative 
functions satisfying 

h'^WPe - (^hPe\\L^h,T;{H''{n)y) + \\Pe\\Li{o,T;m{n)) < C for all h>0, 

where C > is a constant which is independent of e and h. Then (p^) is relatively compact 
m LP^(0, T; U'-{^)) for all i < q and for all r < 2d/{d - 2) %f d > 3, r < oo %f d = 2, and 
r < oo if d = 1. 

A variant of this lemma is due to Dubinskii, see [23, Theoreme 21.1, Chapter 1] for a 
proof. A simple proof is achieved by applying the lemmas of Aubin [27] and Chavent-Jaffre 
[5]. Since the result is of interest by itself, we provide the (short) proof. 

Proof. The function f{x) = x^^^, < x < oo, is Holder continuous with exponent 1/p. 
Therefore, with u = p^, by the lemma of Chavent-Jaffre [5, p. 141], 

IIPelUi/p,2p(f7) = ||/(M)IUi/P,2p(n) < ^hllnf^n) = '^\\P^e\\H?(ny 

This shows that (p^) is bounded in LP''(0,T; W^^P''^p{VI)). By Aubin's lemma [27, Theorem 
6] and the compact embedding W^^P''^p{^1) L'p^{^1) (r is as in the lemma), (pe) is relatively 
compact in Lp\0, T; Lp^'^Q)) for all i < q. D 

The following result, which will be used in this paper, is a consequence of Lemma 6. 

Lemma 7. Let ^7 C M'^ (c? > 1) be a bounded domain with d^l E C°'^ let T > 0, t > 0, 
and let tk = kr, k = 0,...,N, with Nt = T be a decomposition of the interval [0,T]. 
Furthermore, let p > 1, q > 1, and s > 0, and let (pr) be a sequence of nonnegative 
functions, which are piecewise constant in time on (0,T), satisfying 

T'^Wpr - crTpT\\L^T,T;{H''{n)y) + \\p^r\\Li{o,T;m{n)) < C for all r > 0, 
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where C > is a constant which is independent of r. Then (p^) is relatively compact in 
LP^{0,T;LP''{Q)) for all £ < q and for all r < 2d/{d - 2) if d > 3, r < oo if d ^ 2, and 
r < oo if d = 1. 

Proof. Since pr is piecewise constant in time, we can write pr{-, t) — pk for t e ((/c — 1)t, /ct], 
k = 1, . . . ,N, for some functions pk- 

Case h < T. The difference pr — CFhPr partially cancels for h < t, and we obtain, for 
A; = 1, . . . , - 1 and t > h, 

\\Pr[-,t) - {ahPr)[-:t)\\(^H^(^a)y - < Q 

Therefore, by assumption, 

^-1 rtk+h 



\\Pk+i — Pk\\(H''(a)ydt 

||Pfc+l - Pfc||(if^(n))' = / WPri-^t) - {(7TPT){-,t)\\{H'>{Q.)ydt <C 

1 1 1 1 ^ th 



k=l k=l 

uniformly in h < t. 

Case h > T. There exists m eN such that tm < h < tm+i- Then, for t e {tk+m-i, tk+m], 
k^l,...N -m, 



UPt- crhPr){-,t)\\(H^(n)Y 
We compute 



\Pk+m — Pk~l\\{H''{n)y if tk+m-1 < t < tk-l + h 

\Pk+m — Pk\\iH''{n)y if tk-l + h <t < tk+m 



rtm+i ^-"^ rtk-i+h 

/ ||Pm+l — Pl||(i?''(f2))' m Pfe— 1 II 



\Pt — Cr^Pr||Li(/i,r;(i7=(n))') 

\pT — (^hpT\\{H'{Q.)ydt + ^ / IIpt- — ahpT\\{H''{Q,)ydt 

k=2 ^k+m—l 

ydt 

I" I • V--/' 

N—m „+, 

+ / \\Pk+m - Pk\\iH'>in)ydt 

k=2 "^tk-i+h 

N-m 

{tm+i - h)\\pm+i - Pi\\{H^{a)y + {h - tm) ^2 \\Pk+m " Pk-l\\{H'>{a)y 

k=2 

N-m 

+ {tm+1 -h) ^ ||pfc+ 

m Pfc||(i?^(f2))'- 

k=2 
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We employ the estimates h — tm < f, t^+i — h < r and the triangle inequahty: 

m+l 

WPt - crhPT\\L^h,T;{H''{n)y) < ^ X] \\Pj - Pj -iWiH'^in))' 

m N—m m N—m 

+ ''"X^ \\Pk+j - Pk+i-lWiH^in))' + X] \\Pk+j - Pk+j-l\\(H'{n))'- 

j=0 k=2 j=0 k=2 

Since 

N—m m m N—m+j N 

k=2 j=i j=i e=2+j e=2 

for numbers ai > and all < i < m, it follows that 

N-l 

WPt - 0"/iPr||Ll(/i,T;(H=(n))') < 3(m + l)r X \\Pk+l - PkWiH'^iQ))' 

k=l 

= 3{m + l)j \\pr - crrpr\\{H-{n)ydt. 

Thus, using (m + l)r < h + t < 2h, 

II II .. ^, 

WPt — (^hPT\\L^(h,T;(H''(n)y) < ~\\Pt — 0"TPr||Ll(r,r;(//»(n))') < 

We conclude that in both cases, for all h > 0, 

WPt — <^hPT\\L^{h,T;{H''{n)y) < Ch, 

and this estimate is uniform in r > 0. We apply Lemma 6 to conclude the result. □ 

3. Global existence of weak solutions 

In this section, we prove Theorem 1. Let in the following ^7 C M"' (1 < < 3) be a 
bounded domain with dQ G C^'^. The smoothness assumption on the boundary of the 
domain is needed for applying elliptic regularity results. 

3.1. Solution of an approximate problem. We show first the existence of a weak 
solution to an approximate problem which is obtained by semi-discretizing (l)-(2) with 
respect to time and by regularizing the equation for the cell density. For this, let T > 
and K and split the time interval in the subintervals 



K 

(0,T] = [j{{k-l)T,kT], T = T/K. 
k=l 
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For given (pfc_i,Cfc_i), which approximates (p, c) at time r{k — 1), we wish to solve the 
approximate problem in the weak formulation 

1 



{{pk - Pfc-i)0 + a{ck - Ck-i)ip)dx 
+ e {ArkA(j)+\Vrk\^Vrk-V(f) + rk(j))dx = / {pk - Ck)^dx 



n 



where the entropy variables are given by 

and (0,'?/') G H^{Q) x H^{Q) is a test function pair, well defined for n > 1. We prove now 
the existence of a solution to (9) recalling that 

, , m + n — 1 

(10) P= — ^ — . 

Remark 8. For proving the existence of weak solutions, the regularization e J^(ArfcA0 + 
Tkcfydx would be sufficient. The additional term is helpful when deriving energy estimates 
which lead to the uniform boundedness of the solutions to the parabolic-elliptic system, 
see Section 4. □ 

Proposition 9. Let Vt d {1 < d < ?>) he a hounded domain with dQ G C^'^ . Further- 
more, let {vk-i, bk-i) G x with r^-i > m f2 and let m > 0, n > 1 
he such that 1 — 2/d < p < m with p given hy (10). Then there exists a weak solution 
{rk,bk) G H'^{Q) X H^{Q) to (9) satisfying Vk > in Q. 

Proof. Step 1: Formulation of a modified problem. In order to solve (9) in terms of (r, 6), 
we set 

w{r) = p = [ r ) , r+ = max{0,r}. 



n 



We wish to solve first the system 
1 



T 



n 



{{w{r+) - w{rk-i))(j) + OL5{bk - bk-i)ip)dx 



n 



+ e (ArA0 + iVrpVr ■ + r0)dx = / {w{r^) - 6b)^ljdx 



n 



where {(pjip) ^ H'^{VL) x H^{VL). Notice that the assumption p < m is equivalent to 
m — n + 1 > 0, which is needed for the term w{r^)^~^^^ to be well defined. The minimum 
principle shows that any weak solution (r, b) to this problem satisfies r > in f2. Indeed, 
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using (r_, 0), where r_ = min{0, r}, as a test function, and observing that w{r+)'Vr_ = 0, 
we obtain 

-- / w{rk-i)r-dx +— f w(r+)'"""+^|Vr_prfx + £ / ((Ar_)^ + |Vr_|^ + r!.)(ix = 0. 
'T Jn ^ Jo. Jo. 

Since all three integrals on the right-hand side are nonnegative, we conclude that r_ = 
and r > in f2. 

Step 2: The linearized problem. Let a E [0, 1] and (f, 6) G H'^^^{Vl) x be given. 

The Sobolev embedding H'^^^iVt) C^{Vl) for d < 3 shows that w{f^) is bounded. Hence, 
the following linear problem is well defined: 

(12) a((r,6), (0,^')) = for all (</.,7/>) e ^/^(n) x 

where 

+ £ / (ArA0 + |Vf|Vr ■ V0 + r0)(ia; + 5 / bipdx, 
Jn Jn 

f{(j),i)) = -- I {{w{f+) - w{rk^i))(j) + aS(b - bk-i)ij)dx + [ w{f+)i)dx. 



T 



The function a : {H^{Q) x — t- M is bilinear and continuous due to the Sobolev 
embedding H'^{n) C H^/'^{n) C W^'^{n) C ^0(0) for d < 3. Here, we need the assumption 
m — n + 1 > 0. The function / : H'^{Q) x H^{Q) — )■ M is linear and bounded which is a 
consequence of the estimate 

w{rk-i)(pdx < |k(rfc-i)||L"(Q)||0||L"/{"-i)(f7) < C\\(f)\\H^n), 



n 



for some constant C > 0, since r^.i G L"/*^" "'^■'(^) gives w{rk-i) G Moreover, a is 

coercive: 

a{{r,b),{r,b)) = J (^^w(f+)'"-"+^|Vr|2 + 5|V6|')dx 

+ e ((Ar)^ + |Vfp|Vrp + r2)(ia; + 5 / 6^da; 
Jn Jn 

>C{e,S){\\r\\l,^^^ + \\b\\ 

H^{n))^ 

for some constant C{e,6) > 0, since dQ G C^'^ (see Troianiello [28], p. 194). The Lax- 
Milgram lemma now implies the existence and uniqueness of a solution (r, 6) G H'^{Q) x 
H\n) to (12). 

S'tej* 3: The nonlinear problem. The previous step allows us to define the fixed-point 
operator S : [0,1] x H'^/^Q) x L'^^Q) H'^/^Q) x L\n) by S{a,f,b) = (r,6), where 
(r,6) G X iJi(l]) is the unique solution to (12). It holds S{0,f,b) = (0,0) for all 

(f , 6) G //^/^(fi) X L^(f2). Standard arguments prove that S is continuous and compact, 
taken into account the compact embedding of H'^{n) x H^(n) into H'^/^iVt) x L'^iVt). 
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It remains to show that there exists a constant C > such that for all fixed points 
(r, 6) G H'^/'^iyt) X and cr G [0, 1] satisfying S{(r,r,b) = {r,b), the estimate 

(13) \\{r,b)\\Hr/i(n)xL'Hn) < C 

holds. Let (r, 6) be such a fixed point. Let us first assume that cr = 1. Then (r, 6) is a 
solution to (9). By the first step of the proof, we have r > in f2. Moreover, we can easily 
derive a uniform bound for p = w{r) by employing (1, 0) as a test function in (9): 



pdx = / Pk^idx — TE rdx < / pk-idx, 
1 Jq Jn Jq 

since r is nonnegative. By iteration, we infer that 
(14) IIpIUhh) < IIPollLi(n)- 

The uniform estimate (13) is a consequence of the following discrete entropy estimate, 
which settles the case a = 1. The case a < 1 can be treated similarly. 

Lemma 10. Let {r,b) G H'^iVl) x H^iVl) be a solution to (9) and let 1 — 2/d < p < m. 
Then 

rmn ..^ „,,2 n 112 

+ £^{\\M\hin) + l|Vr||i4(f^) + ||r||i2(f^)) < S(pfc_i, Cfe_i), 

where the entropy E(p,c) is defined in (6), p = w{r) = ((n — l)r/nY^^'^^^\ c = 5b, and p 
is defined in (10). 

In order to prove this lemma, we employ the test function (r, 6) = {np^~^/{n — l),c/6) 
in (9): 

^ ^ -p"~^{p- Pk-i) + ^c(c- Ck-i))dx 



r /o — 1 6 



(15) + / C^lVp^l^ + -J\Vc\^ + c^))dx + e [ {{Ar f + \Vr\'' + r^)dx=\ I pcdx, 

where p^-i = ui{rk_i) and c^-i = Sbk-i- Since n > 1, the mapping g{x) = x^, a; > 0, is 
convex, which implies the inequality g{x) — g{y) < g'{x){x — y) for all x, y > 0. Hence, the 
first integral on the left-hand side of (15) is bounded from below by 

- / (^ip" -pU) + i^Ac'-cU))dx = -iEip,c)-Eip,^^,c,^,)). 

For the estimate of the right-hand side of (15), we employ first the Holder and Young 
inequalities: 



1 /■ . - 1 1 ,2 1 



- / pcdx < T||p||L«(f7)||c||i,'(n) < ^||p||M(n) + ^\\4m{n), 
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where q > 6/5 if d = 3, 1 < g < oo if d < 2, and q' = q/{q ~ 1). In the last step, we have 
used the continuous embedding H^{Q) ^ L'' (Q) which is vahd since g' < 6 if (i = 3 and 
g' < oo if < 2. By Lemma 5, we find that 

1 f TflTL 1 

- j^pcdx < ^||Vp^||i2(f^) + ||p||Li(n)) + ^||c||^i(f^). 

The assumptions of the lemma are clearly satisfied for d < 2. If d = 3 we can choose 
q = 2d/{d + 2) = 6/5 for p < 1 and q = p + 1/2 > 6/5 for p > 1. Putting together 
the above estimates and the bound (14), this finishes the proof of Lemma 10 and of 
Proposition 9. □ 

3.2. Uniform estimates. Let (r^, b^) be a solution to the approximated problem (9) and 
set pk = w{rk), Ck = Sbk- We define the piecewise constant functions 

(p(^'^),r(^'^),c(^'^))(x,t) = {pk,n,Ck){x) foTxen, te {{k - l)r, fcr]. 

We denote by D^p{-,t) = (p(-, t) — p{-,t — r))/T the discrete time derivative of p(-, t), where 
t > T. In terms of the variables (p*^'^'^\ c'-^'^-'), system (11) can be formulated as 

0= / {DrP^''^\(^)dt+ [ [ (V(p("'"))"-p("'")Vc("'")) ■ V0(ia;(it 
Jt Jt Jq 

(16) +e [ [ (Ar(^'^) A0 + | Vr^^'^^ Vr^"'^) ■ V0 + r^^'^V)^^^ dt, 

Jt Jn 

= a [ {D^c^''^\ij)dt+ I I (Vc("'") + 5V(p("'"))") ■ V^dxdt 

Jt Jt Jq 

(17) + / / (p("'") -c("'"))V^c/xc/t 



for all smooth test functions and ip, where (■, ■) is a dual product. We set Qt = fix {0, T) 
for given T > 0. Before we can perform the limit (e, r) — )• 0, we need to prove some uniform 
bounds in e and r. The following result is a consequence of the discrete entropy estimate 
of Lemma 10 and the bound (14), after integrating with respect to time. 

Lemma 11. Let T > and 1 — 2/d < p < m. Then the following uniform bounds hold: 

(18) l|P^'^'^''l|L°°(0,T;Li(n)nL"{Q)) + II V (p^^'^"* )^ || l2 (Qj,) < C, 

(19) V^\\r^''^^\\mo,T;mm + ^||Vr(^'^)|U4(f,^) < c, 

(20) tt||c^^'^^ 11^00(0,^.^2(5^)) + ||c^^'^''||L2(o,T;Hi(n)) < C, 

where C > is here and in the following a generic constant independent of e and t, and 
p is defined in (10). 

Under additional assumptions on the exponents n and p, we are able to derive more a 
priori estimates. 
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Lemma 12. Let p < min{m, n} and Q := n / d + p > 1 and set 

^1 = 7^ e 1,2, S2 = y—^ G 1,2, ss = -^e{l,2). 

Q+m-p Q+n-p Q+l 

Then the following uniform bounds hold: 

(21) \\{p^'''^r\\mo,T;mm + ||p(^'^)|U.Q(n,) < C, 

(22) |l(P^^'^'')''"llL''i(0,T;iyi.»i(n)) + ll(p''^'^'')"llL''2(0,T;Tyl'''2(Q)) < C, 

(23) ||p(^'")Vc(^'^)|U^3(n,)<C, 

(24) ||-DrP^^'^''||LS(T,T;{/f3(Q))/) + a || L)^C^^'^^ || Ls(^,T;(H3(n))') < C, 

where s = min{si, S2, S3}, s = min{s,4/3}, and p is defined in (10). 

We remark that the condition Q > 1 is equivalent to p > 1 — n/d, which in particular 
implies the condition p > 1 — 2/d a.s explained in the introduction after (7). 

Proof. We set p = p^^''^^ and c = c*-^''^-' to simplify the notation. By the Poincare inequality, 
we find that 

mhinr)= fmhmdt<C f \\Vp^\\l.^^^dt + C fmlr^n)dt. 
Jo Jo Jo 

Since p is uniformly bounded in L°°{0,T] L"(f2)) and since we have assumed that p < n, 
the right-hand side of the above inequality is uniformly bounded. This shows the uniform 
bound for p^ in L^(0, T; H^{fl)). Next, the Gagliardo-Nirenberg inequality with 9 = p/Q < 
1 gives 

||.||2Q _ ||^||2Q/P . ^ /"'^ II .||2Qe/p|| p||2Q(l-e)/p , 

•/ 



lL°°{0,T;L"(n)) / llA^ WH^U) 
J 

This shows (21). 

For the proof of (22), we observe that Si > 1 is equivalent to n/rf + n > 1, which is true 
since n > 1, and that Si < 2 is equivalent to p < m, which holds by assumption. Hence 
Si G (1,2]. Let first Si < 2. We apply the Holder inequality with exponents 7 = 2/si>l 
and y = 2/(2 - Si): 

I|V(P")|I1^..(^,) = i^T r [ p^-~'^'^\Vp^rdxdt 

^ Jo Jn 

because of (21). If si = 2, it follows that m = p, and the conclusion still holds. The 
estimate for p"^ is shown in a similar way by applying the Holder inequality with exponents 
7 = 2/si > 1 (if S2 < 2) and 7' = 2/(2 - si) to p"" = p'^'^ff. Hence, p™ is uniformly 
bounded in L'^{0,T; W^'''{n)). 
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The estimates for V(p") and p" in L^'^{Qx) are proved analogously by replacing si by 
52- The relation S2 < 2 is equivalent to p < n, and S2 > 1 is equivalent to n/d + m > 1, 
which is true since n/d + m>n/d + p> Iby assumption. This proves (22). 

Estimate (23) is again a consequence of the Holder inequality, with exponents a = 2/s3 
and a' = 2/(2 - S3): 

||pVc||L»3(nj,) < ||p||L2s3/(2-.3)(f^^)|| Vc||l2(Qj,) = ||p||i2Q(f^y)|| Vc||i2(j^y) < C. 

We turn now to the estimates for the discretized time derivatives. Let G L*'(0,T; 
H^{il)), where s' = s/{s - 1) > 2 and s' = max{s',4}. Then, using 1 < s < 2 and (19), 
(22), and (23), 

-T 

{Dt-p, (j))dt 



I I {Vip"")- pVc) ■V(t)dxdt-e f [ (ArA0 + | Vr|^Vr ■ V0 
Jt Jn Jt Jn 



r(j))dxdt 



< (l|V(p^ 



+ ||pVc|U.(n,))||V0|| 



+ 4M\L\nT)\\^(P\\L\nT) + ^l|Vr||i4m ^||V0|U4(^^^) + e\\r\\L2^nT)U\\L^ 



^ C\\<P\\LS'{0,T;H'-^{n))- 

Furthermore, since p G L'^^{Qt) C L'^IQt), we obtain 



{D^-c, (f))dt 



T 



< 



(Vc + 5V(p")) • V<pdxdt + j I (p- c)(t)dxdt 
(||VclU.(n,) + 5||V(p")|U.(n,))||V0|L.'(^^) 



+ (||pI|l2{Qt) + l|c||L2(nT))ll0lU2(nT) 
< C'II0|Il'''(O,T;H3{Q))- 

using (20), (21), and (22). Thus we have proved (24). 



□ 



3.3. The limit of vanishing approximation parameters. We show first the strong 
convergence of the sequence (p*-^'^^). 

Lemma 13. Let the assumptions of Lemma 12 hold. Then, up to a subsequence, 
(25) p(^'^) p strongly m ^^'^(0, T; ^^^(fi)), 

for all r < Q, where Q = n/d + p > 1. 

Proof. First, consider p < 1. We show that (p*-*^'"^^) is bounded in L^(0,T; W^'^{^1)), where 
£ = 2Q/{Q + 1 — p) > 1. Notice that 1 < i < 2 since the former inequality is equivalent to 
Q > 1 — p which is true since Q > 1 by assumption, and the latter property is equivalent 
to p < 1. We apply the Holder inequality to Vp*^^'^-* = (l/p)(p'^^'^))^~^V(p*^^''^^)^ yielding 



|Vp(^'-)|U.(n,) < c||(p(^'^))i-^|L2./(2-.,(n.)l|v(p(^'^)^^i 



lL2*/{2-<)(nj,)| 

I 
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since 2£(1 — p)/(2 — £) = 2Q. Thus, taking into account the bound for (p^^'"^^) in L'^^iVlT)^ 
by (21), this proves the desired bound. 

Next, we claim that the embedding W^'^{VL) =— )■ L^*5(f2) is compact. This is the case if 
1 — d/l > —d/{2Q), which is equivalent to 2n/d + m > 1. In order to show this inequality, 
we observe that the assumption p > 1 — 2/d is equivalent ton + m>3 — A/d. U d = 2, 
this implies that 2n/d + m = n + m > 1; if d = 3, we find that 

2n 2, ^ m 2, ^ 2 5 

d 3^ ' 3 3^ ' 3 3 

Hence, 1 — d/i > —d/ {2Q) is satisfied for d < 3. In view of the bound for the discrete time 
derivative of (p*-^'^^), see (24), Lemma 7 (take p = 1 in the lemma) implies the existence of 
a subsequence of {p^^''^^) (not relabeled) such that (25) holds. 

Finally, if p > 1, we can apply Lemma 7 to conclude the strong convergence result. □ 

The estimates of Lemmas 11 and 12 imply the existence of a subsequence of (p*-^'^-*), 
which is not relabeled, such that, as {e,T) — )■ 0, 





^ C 


weakly 


in 


L\0,T;H\Q)), 




Zl 


weakly 


in 






Z2 


weakly 


in 






Z3 


weakly 


in 




DrP^''^^ 


dtp 


weakly 


in 






■ dtc 


weakly 


in 


L'{0,T;{H\Q)y) if a 







weakly 


in 


L\0,T;H\Q)), 




^ 


weakly 


in 





The limits of the nonlinearities are easily identified since (a subsequence of) (p*^^'"^^) con- 
verges strongly in L^'5(f2j'), where Q > 1. Hence, up to a subsequence, p^"^'^) — )■ p a.e. in 
VLt and (p^^''^))"^ — > p™, (p^^'"^)) — )• p" a.e., implying that Zi = p™, Z2 = p". Moreover, the 
strong convergence of (p*^^'"^^) in ^^(^2^) and the weak convergence of Vc'-'^'^-' in L^(fiT') give 
the weak convergence of (p*^^''^^Vc^^''^'') to pVc in L^{ilT), implying that = pVc. 

The above convergence results are sufficient to pass to the limit (e, r) — )■ in (16)-(17) 
leading to (l)-(2). The Neumann boundary conditions are satisfied in the weak sense, and 
the initial conditions hold in the sense of L*(0,T; (if^(il))'). Since in the limiting equation 
the regularizing terms vanish, test functions in (0,T; W^''^ {VL)) are sufficient to obtain 
the boundedness of the diffusion and drift terms. A density argument now completes the 
proof. 

4. The parabolic-elliptic system 
The parabolic-elliptic system corresponding to (2) is given by 
(26) dtp = Ap™ - div(pVc) , = Ac + 5Ap" + p - c. 
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subject to the no-flux boundary conditions 

(27) (Vp™ - pVc) ■u = 0, V(c + (5p") ■ z/ = on dn, t > 0, 
and the initial condition 

(28) p(-,0)=po in a 

Similarly as in [14] we introduce a new unknown corresponding to the diffusion terms in 
the second equation, v = c + Sp"', and rewrite system (26) in terms of p and v. 

dtp = Afp^ + ^^^p^+i) -div(pVt;), 

(29) V n + V J ^'^ ^' 

= Av + p + dp"" -V, 

subject to the no-flux boundary conditions and initial data (28). In the case n = 1 one 
simply obtains the Keller-Segel model with nonlinear diffusion, which is known to prevent 
blow up [4, 21, 22]. Here, the situation is different, since for n > 1 we obtain additionally 
a nonlinear growth term in the equation for v. We will show that nevertheless this system 
satisfies the properties used by Kowalczyk in [21] to obtain an L°°{Qt) bound for p. The 
difference here is that instead of the conservation of mass we have to make use of the 
uniform boundedness of ||p||Ln(n), which only holds for finite times. 

The main advantage of the parabolic-elliptic system is that it allows for another entropy, 
since we will show that we can use powers of p as test functions in the elliptic equation, 
see Remark 2. 

Lemma 14. Let {p,v) be the solution of the parabolic- elliptic system (26) constructed in 
Theorem 1 with po G L"{VL) and po > 0. Then 

p > 0, f > a.e. in VLt, 

and p satisfies the additional entropy estimate 

(30) sup /^c^x+^l|V/||i.(^^) + |||Vp"||i.(,^)<C. 
tG{0,T) Jn ^ ~ ^ P 

Proof. We start again from the regularized problem to derive the nonnegativity of v and 
the additional entropy estimate rigorously. As in Section 3 we skip the index k: 

[ Drp(pdx = - [ vfp™ + 5^^p"+^V [ pVv-V(t)dx 

Jn Jn ^ n + 1 / 

-e {ArA(p+\Vr\^Vr -Vcp + rcfjdx, 
Jn 

= - / Vv ■ Vipdx + / (p + Sp^)iidx - I vipdx, 
Jn Jn Jn 

for appropriate test functions and ip. The existence of a global weak solution is proven 
in Proposition 9. The a priori nonnegativity of v follows from a standard argument by 
testing the equation for v with v~ = min{0,f}, see, e.g., [21]. The nonnegativity is clearly 
preserved when performing the limit of vanishing parameters. 



(31) 



n 
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To derive the additional energy estimate we first use r = ^^p""^ as a test function in 
tlie equation for p, integrate in time, and insert the elhptic equation for v: 

< / V't>.Vp"di*-£(||Ar||i.,„^, + ||Vr||l.,„^, + ||r||i,,„^l) 

< 1 Vv-Vp'^dxdt= I p^'+^dxdt + S I p^'^dxdt-f vp^'dxdt 
<C^+^l|Vp1li.(n,). 

Here, we have used the nonnegativity of p and v and the uniform boundedness of ||p||L"(n) 
(with respect to k) together with the Gaghardo-Nirenberg inequahties in the following way: 

— II <r r'iiY7n"l|("+^)^i/'"ll /i"l|("+^)(^^^i)/" _L II /i"l|("+^)/" 
P aj- — IIP ||^(„+i)/„(j^) ^ II vp 11^2(5^) IIP 11^1(0) "t" IIP llLi(n) 

< C^I|Vp"1li"|5^^/" + ^ ^ ^l|Vp"||i.(n) + C{5\ 
where di = 2d/{{n + 1)(2 + d)) G (0, 1), and for 62 = d/{d + 2), we obtain the bound 

[ p'-dx = Wp-Wh^n) < C||Vp1lf4^)||p1li^f^;^^ + mlnn) < ^l|Vp"||i.(f,) + 

Since the constant C{6) is independent of (e, t), we can perform the limit of vanishing 
parameters, which completes the proof. □ 

Theorem 15. Let po G L'^{Q), where Q G M.'^ is a bounded domain {1 < d < 3) with 
dQ G C^'^, and let the assumptions (5) hold. Then the parabolic- elliptic system (26)-(28) 
has a global weak solution satisyfing 

||p||L°°(0,T;L°°(n)) + ||c||L°°(0,T;L°°(n)) < C (T) 

for any T > 0. 

In particular, test functions in L^(0,T; H^{Q)) are admissible for (26). 

Proof. The iterative method of Alikakos used by Kowalczyk to derive an L°°{Qt) bound 
for p requires test functions of the form p'' for g > n as well as the uniform boundedness of 
Vv. Thus, we prove first that Vv G L°°{Qt) and p G L°°(0,T; L«(fi)) for suitable g > 1. 

Step 1: Proof of V G L°°(0, T; ^^l'°^(^])). If p G L°°(0, T; L3"+i(fi)), then elliptic regu- 
larity for 

—Aw + v = p + 6p^ in Vt, Vf ■ z/ = m. dVt 
implies that v G L°°(0, T; VT^.s+i/nj^f^^^ (^^^^^ g_g_^ ^^^j^ p_ ^26), and hence, by Sobolev 
embedding, v G L°^(0, T; Vri'°^(fi)). In order to show that p G L°°(0, T; L3"+1(J1)), we 
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employ in the regularized equation (31) the test function (3?t, + l)p'^'^ and integrate in time: 



[p'"+'dx~ [ p3"+Mx + 2a||Vp(™+=^")/li.(^^) + 25611 Vp2"+^ 
Jn Jn 

+ (3n + l)e (^j ArAp3" + | Vr^Vr ■ Vp^" + rp^"^ dx dt 

<3n Vp^"+^ • Vvdxdt 
Jnr 

= 3n [ p^^'+^dxdt + 35n I p^'^+^dx dt - 3n I p^'^^^vdx dt 

Jut J^t J^t 



/2||2 



< A5n / p'^'^+'dxdt + C, 

where we have used the nonnegativity of p and v, Holder's inequality, and 3n + 2 < 4n + 1. 
Moreover, we have set a = 6mn{3n + l)/(m + 3n)^ and b = 6n^(3n + l)/(4n + 1)^. We 
apply the Gagliardo-Nirenberg inequality with 9 = d/{d + 2) and Holder's and Young's 
inequalities to estimate 



/ n^^+l^-r - Iln2"+1/2||2 ^ ^ || 2n+l/2 || 25 || 2n+l/2 || 2(1-9) n 2n+l/2||2 

J^P ax — \\p ||i2(f^) ^ u II vp \\l^{q)\\P iiL^n) WP Wl^ 
The last summand is estimated by interpolating, for any /3 > 0, 

2n+l/2||2 



IP 



Ihiu) < IIP^'""''^/'lli.(n)IIPlUMn) <C [ p'^'^'dx <P f p'-^'dx + 



since p G L°°{0,T] L^{Q)). Finally, combining the above estimates and by choosing /3 
appropriately, we obtain 

/ p'-^'dxdt<Uvp'-^'/YLHn,) + C. 
It remains to bound from below the terms arising from the regularization: 
e [ ( Ar Ap3" + | Vr | ^ Vr ■ Vp^" + rp^") dxdt 

>3e f (p^''^UAr? + 2(l + —]p''+^\Vr\^Ar + p^^+'\VrAdxdt 



3e j (^p2"+i(Ar)2 + 2(^l + ^)p"+2|Vr|2Ar + p2"+i|Vr|4^ dx. 
3e j (^p"+i/2Ar + (l + ^p^l^\Vrij dxdt 



+ 3e I p' (^p2("-i) - (^1 + -Lj' ) |Vr|^dxrft 
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>e [ p'( - (l + i-)''] \Vr\^dxdt 



{p"-i<l+l/2n} 



> -eC\\Vr\\l,,^. > -C . 



where C is independent of {e, r) due to the basic entropy estimate for the regularized 
system. Thus, we obtain in the hmit (e, r) — )■ for the weak solution of (26), 

Jn 

This implies that p G L°°(0, T; L3"+^(fi)), which proves the claim. 

Step 2: Test functions in L'^{0,T] H^{Q)) are admissible. We have to verify that Vp™" G 
L^(0, T; L^(n)). We recall that the restrictions on the exponents are n — 1 < m < n + 1 or, 
equivalently, p<'m<n + l<2n + l/2. Therefore, in view of Step 1, we can interpolate 



p2 



I \Vp^\^dxdt= [ p'^^"'-P^\VpP\'^dxdt 

< [ \VpP\^dxdt+ I p2(2"+l/2-p)|VpP|2rfa;rft 

J{p<l} J{p>l} 

<l|v/||i.(^,) + c7||vp^"+^/^lli.(^^)<a 

Step 3: We now proceed to make the estimates derived by Kowalczyk rigorous. To this 
end, we will use powers of the cut-off functions px = min{p, K} as test functions. Due to 
(30), p^ G L^(0,T; H^{Q)) for any q > n and according to Step 2, it is an admissible test 
function. Let us introduce the notation 

w{x, t) = {pk — k)^ for some A; > 0. 

We test the equation for the cell density with (g + l)w'^, where q > n: 

dt [ w'^-^^dx = -5{q + l) I Vp"" -VwHx-iq+l) [ Vp'^-Vw'^dx 
Jn Jn Jn 

+ {q + l) / pVw" ■ Vvdx 
Jn 

< -5n(g + 1)A;"-^ [ V p ■Vw'^dx - mq{q + I) [ p'^-^w''-^\V p\^dx 
Jn Jn 

+ 2g||Vt;|Uoo(n) (^J w^'^+'^/^\Vw^''+'^/^\dx + k j w^^-^^/^\Vw^'^'-^^'^\dx^ 

Neglecting the second term on the right-hand side and employing Young's inequality to 
the last two terms, we arrive at 

dt f w^^'dx <-6^^k''-' I \Vw^''^^^'^\^dx + Cq{q+l) I w'^^^ dx + C qiq + I) , 
Jn (1 + ^ Jn Jn 
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where C > depends on, e.g., 6 and || Vf ||loo(q). Starting from this inequahty, Kowalczyk 
[21] employed the iterative method of Ahkakos to obtain 

\\w{-,t)\\L--{n) < C{\\wo\\Lo^(n)) < C(||po||L->(n)) for all < t < T. 

Hence, for po ^ L°°{^t), we have w = {px — ^)+ G L°°{Qt) foi^ any K, k > with an 
independent bound for the norm. We can let — ?■ oo to deduce 

||p(-,t)||L-(n) < C(||polU-(n)) for all < t < T, 
which finishes the proof. □ 



5. Long-time behaviour 

The (modified) Keller-Segel system possesses the constant homogeneous steady state 
p* = c* = M/meas(f2). Let us consider the following system, which is equivalent to 
(l)-(2): 

(32) dtp = Ap"' -divipVc), 

(33) adtc = Ac + 5Ap" + p- p* -(c-c*) in n, t > 0. 

In the case of linear diffusion terms, m = n = 1, the decay of the relative entropy yields 
the convergence of the solution towards the homogeneous steady state for large enough 6 
if d = 2 [14]. The corresponding relative entropy for the system under consideration with 
nonlinear diffusion is 



E*{t) 



{p-p*Y {c-c 



\2 



+ a- 



n-l 25 



dx. 



Notice that the nonnegativity of E* is only guaranteed if n is an even integer. In particular, 
E* is not well defined for general real n > 1, since p — p* may be negative. Formally, testing 
(32)-(33) with (n(p — p*Y~^ /{n — 1), (c — c*)/5) we obtain the evolution equation for the 
relative entropy 

^ + mnj^p--\p - pT-'\^P?dx + \\\c- c*||^.(^) 

= n f p{p - p*Y~^V c ■ y pdx - n [ p'^-^V p ■ Vcdx + ^ [ {p - p*){c - c*)dx. 
Jq Jn Jq 

We see that for general values n > 1, also the entropy dissipation terms are not necessarily 
nonnegative. Moreover, the chemotactic drift term and the term arising from the cross- 
diffusion perturbation cancel out only if n = 2. These comments motivate us to consider 
the case n = 2 only. We recall that, due to (5), n = 2 implies that 1 < m < 3. Thus we 
restrict us to the special case 

m = 1, n = 2, 

for which the global existence of a weak solution is guaranteed. This choice allows us to 
show the decay of the solution to the homogeneous steady state for certain values of 6. For 
the fully parabolic system we need additionally a smoothness assumption on the solution, 
since the weak solution, obtained in Theorem 1, cannot be used as a test function (in 
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contrast to the parabolic-elliptic system), see Remark 2. Notice that we cannot start from 
the regularized problem, since there the mass of p is not conserved, hence the system does 
not possesses a constant homogeneous steady state. 

Now, we are in the position to prove Proposition 4 (see the introduction). 

Proof. Step 1: Decay of the relative entropy. We wish to employ (2(p — p*), (c — c*)/25) as 
test function in (32)- (33), which is allowed if a = 0. On the other hand, when a = 1, we 
need the smoothness assumption. 

Since n = 2, the evolution equation for the relative entropy 

{p-pr + a- 



E*{t) 
reduces to 
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dx 



The difference to the entropy estimate in the existence proof is that we can now apply the 
Poincare inequality to p — p* to derive the decay of the relative entropy for certain values 
of (5 : 

^ jj.p- Plic- c*)dx < ^\\p - p*\\L2in)\\c - c*\\L2(n) < ^||p - P*||i2(n) + ^IW - c*\\l2^^) 

<^l|Vp||i.(^) + ^||c-c*||i.(^). 
Finally, we obtain the entropy estimate 

j^E*{t) + (2 - ^) II Vplli.(^) + ^||c - c*||l,.(^) < 0. 
We set K = min{l, 46 - Cj,}/{A6) > 0. Then 

—E* < -2kE*, 
dt 

implying the exponential decay of the relative entropy 

E*{t) < E*{0)e-^^\ t>0. 

If a = 1 this immediately implies the desired decay of p — p* and c — c* in L^(^7). For 
a = 0, the relative entropy only gives the decay of p — p* in L^{Q). 

Step 2: Decay of c — c* for a = 0. Setting v* = c* + 6{p*y, we find that 

y-y* =c-c* + 6{p^ - {p*f) =€-€* + 26p*{p - p*) + 6{p - p*f. 

Replacing c — c* in the elliptic equation for f , 

Q = + p- c = /\v + {p- p*) - {c- c*), 

it follows that 

(34) = A(t; - V*) + (1 + 26p*){p - p*) + 6{p - p*f - {v - v*). 
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Hence, 



Jn Jn 



To determine the decay of ||f — f * ||_f/i(n) we shall derive a uniform bound for v—v*. We prove 
in Step 3 below the boundedness of p — p* in L°°(0, oo; L'^{Q)). Then elliptic regularity for 
(34) gives w-u* e L°°{0,oo]H^{Q)) and, by Sobolev embedding, 'y-'?;* e L~(0, oo; L°°(fi)). 
Hence, proceeding with the above estimate and using Young's inequality, we obtain 

11^ - v*\\min) < ^11^ - + C\\p - P*||i2(f^) + S\\v - t;*||L-(n)||p - P*|li2(o), 

such that 

lit; - v*\\l.^^^ < 2{C + 5\\v- v^h^t^n)) \\p - P*\\h^n) < Ce-^^\ 

This implies for the original unknown c — c* that 

l|c - c*|Ui(o) < C{\\v - v*\\l2(^u) + 2(5p*||p - p*||l2(q) + Hp - P*\\l2(n)) < Ce-'^K 

Step 3: Proof of p - p* G L°°{0, oo; L^(fi)). Using 4(p - p*Y as a test function in the 
first equation of (29) and employing the second equation in (29), we infer that 

dt [ {p-p*ydx + 12 [ {p-p*f\Vp\^dx + 86 [ p\p - p*f\Vp\^dx 

= 12 p{p - p*fVp ■ Vvdx 
Jn 

= 3 [ V{p- p*Y -Vvdx + Ap* [ V{p~ p*f -Vvdx 
Jn Jn 

= 4(1 + 25p*)p* /" (p - p*Ydx + (3 + lOSp*) f {p- p*fdx + 35 [ (p- p*fdx 
Jn Jn Jn 

{v - v*){p - p*Ydx - 4p* ! {v- v*){p - p*fdx. 



The last two terms are bounded by using the nonnegativity of p and v and the Cauchy 
Schwarz inequality: 



{v - v*){p - p*Ydx - 4p* / {v - v*)p{p - p*fdx + A{p*f / {v - v*){p - p*Ydx 

Jn Jn 

<3v* [ {p - p*Ydx + Ap*v* [ p{p- p*fdx + 2{p*f [ {v-v*fdx 
Jn Jn Jn 



+2{p*r {p-prdx. 

Jn 

Together with the estimate resulting from (34), 



1^ - v^WW) < 2(1 + 25p*)' /(P - Pl'dx + 26' / (p - p*rdx, 

Jn Jn 
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we obtain by interpolation: 

dt f {p-p*Ydx + 12 [ {p- p*y\Vp\^dx + 86 [ p\p - p*f\Vp\^dx 
Jn Jn Jn 

<c [ ((p - p*r + (p - p*r + (p - p*f) dx+36 [{p- p*fdx 

Jn Jn 

<C{S) [ ip~p*Ydx + 4:S [ {p-p*fdx. 
Jn Jn 

We already know the decay of ||p — p*\\l'^(q,)'^ hence, it remains to bound ||p — p*||l6(q) 
in terms of the entropy dissipation. To this aim, we use the Gagliardo-Nirenberg with 
6 = d/{d + 2) and the Young inequality: 

IKp - PlYmn) < C{\\Vip - p*)ir.(^)||(p - P*?C[n) + IKp - P*niin)) 
< C\\p{p - p*)Vp\\h^n) + C\\{p - P*)Vp||i.(^) + ^ (i " ''*''^'^) 

<C||p(p-p*)Vp||i.(^) + 6||(p-p*)Vp||i.(^) + C / {p-p*rdx [ {p-p^dx. 

Jn Jn 

The decay ||p — p*\\L'^(n) < Ce~'^^ implies that 

dt [ {p-p*Ydx + Q [ {p- p*y\Vp\^dx + 46 [ p\p - p*Y\V p\^dx 
Jn Jn Jn 

< + / (p - p*Ydx, 

Jn 

and we conclude that p — p* E L°°(0, oo; L'^^Q)), which completes the proof. □ 

6. Numerical Simulations 

This section is intended to illustrate numerically the solutions to the fully parabolic 
system in two and three space dimensions. We compare the results obtained for 5 = and 
6 = 0.005 with various values for the exponent n in (2). The simulations were carried out 
using the COMSOL Multiphysics package with quadratic finite elements. The numerical 
solutions are for illustration only; a more detailed comparison is the subject of future work. 
We choose Q = -Bi(O) for simplicity. 

The two-dimensional case. We consider the fast-diffusion case m = ^ and prescribe 
the initial data 

Po{x,y)=80{x^ + y^-l)\x-0.lY + 5, co{x,y) = for {x,y) E B^{0) 

with M = f^po{x,y)dx = 257r/3 > 8n (see Figure 2 left). The maximal density is 
Pmax = 21.5. We recall that solutions to the classical parabolic-elliptic Keller-Segel model 
(m = 1 and 6 = 0) blow up in finite time when the initial mass M is sufficiently large. 
More precisely, in the radial case, under an additional assumption on the second moment, 
the solution blows up if M > Svr [25] or, in the non-radial case, ii M > An [24]. Since 
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m > 1 leads to global existence results for the parabolic-elliptic system [4, 21], one may 
conjecture that the cell density of the parabolic-parabolic model blows up in finite time 
for M > Stt if m < 1. We confirm this conjecture numerically for the case m = ^ and the 
above initial datum. 




(a) t = 0, p,„ax = 21.5. (b) t « 0.15, p,nax ~ 1212. 



Figure 2. Cell density with 6 = 0. 

The nonlinear diffusion terms cause numerical difficulties whenever the solution becomes 
close to zero. Indeed, the numerical solution may become negative, and the simulations 
break down. Clearly, this can be handled by developing a positivity-preserving numerical 
scheme, similarly as for the porous-medium equation. Since we are using the black-box 
solver COMSOL Multiphysics, we solve this problem simply by a projection method, i.e., 
we replace diffusion terms by A(max{p, e:})™ with m = ^, e = 0.005 and A(max{p, 0})", 
respectively. 

The cell density of the Keller-Segel model with m = ^ and 5 = at time t ~ 0.15 is 
depicted in Figure 2 (right). Shortly after that time, the simulations break down which 
indicates the blow up of solutions. Surprisingly, the singularity forms in the interior of the 
domain in contrast to the classical Keller-Segel model (m = 1 and 6 = 0) for which blow 
up occurs at the boundary. Our numerical experiments confirm this behavior for the model 
with m = 1 and 6 = (results are not shown). Thus, the unexpected behavior seems to 
be an effect of the fast cell diffusion. 

Next, we turn to the case m = ^ and 6 > 0. Figure 3 shows the cell density at time 
t = 1000 for various exponents n = |,y,|,2. The solutions have essentially reached 
their steady state at t = 1000. Notice that, according to (5), the admissible parameter 
range for n is | < n < |. Although some of the values for n used in the simulations 
are theoretically not admissible, the solution exists numerically for all time. However, we 
observed numerical difficulties for large values for n (e.g. n = 6) which may indicate that 
the upper bound for n in terms of m is more than just a technical assumption. We see 
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Figure 3. Cell density at time t = 1000 with 6 = 0.005. 

In the limit of vanishing additional cross-diffusion 5 — )■ 0, we expect that the solutions 
converge to the solution to the corresponding Keller-Segel model with 6 = 0. This is 
numerically confirmed in Figure 4. For 6 = 0.005, the cell density reaches its maximum 
at the boundary (see Figure 3), whereas the maximum is attained in the interior of the 
domain for very small values of 6. Thus, it seems that the cross- diffusion regularization 
produces a stationary state which is more concentrated on the boundary. 

The three-dimensional case. We consider the linear case m = 1 and prescribe the 
initial data 



Po{x, y,z) = 10 + 80{x^ + y^ + z^-l)\x-0Ay, co{x,y,z) = 0, {x,y) E B^{0), 
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(a) S = 10-\ pn,ax « 107. (b) S = 10-6, p„,ax ~ 681. 




(C) 6 = 10-8, p,„ax « 1201. (d) S^O, p„,ax ~ 1212. 

Figure 4. Cell density at time t ^ 0.15 with n = 3/2. 



see Figure 5 (left). In three space dimenions, even for the parabolic-elliptic system, there 
is no critical threshold known for the occurence of blow up. A more complicated functional 
relation between the second moment and the L^/^ norm of the cell density has been derived 
in [7] as a sufficient condition for finite-time blow up. We observed that the numerical 
solution to the fully parabolic model with 5 = breaks down after time t ^ 0.46, which 
may indicate a blow-up formation (see Figure 5 right). 

The initial data in Figure 5 is represented using slices, since the highest values occur 
inside the domain, whereas for the following simulations, we use the level-set representation 
which is more appropriate for demonstrating aggregation phenomena on the boundary. 

In Figure 6 we compare the cell density at time t = 1000 with 6 = 0.005 and n = 
|, j|, |, |. At t = 1000, the solutions have essentially reached the (non-homogeneous) 
steady state. As we already proved in Proposition 4, when performing the same simulation 
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(a) t = 0, pniax « 46. (b) t « 0.46, p,nax « 5460. 



Figure 5. Cell density with 5 = 0. 

with n = 2, the solution converges to the homogeneous steady state. As in the two- 
dimensional situation, the maximal cell density is achieved at the boundary. 
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